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The Cusp/Core problem: supernovae feedback versus the 
baryonic clumps and dynamical friction model 

A. Del Popolo 1,2,3 • F. Pace * 1 2 * 4 * 


Abstract In the present paper, we compare the predictions 
of two well known mechanisms considered able to solve 
the cusp/core problem (a. supernova feedback; b. bary¬ 
onic clumps-DM interaction) by comparing their theoret¬ 
ical predictions to recent observations of the inner slopes 
of galaxies with masses ranging from dSphs to normal 
spirals. We compare the a-V m t and the cr-M* relation¬ 
ships, predicted by the two models with high resolution data 
coming from (Adams et al. 2014; Simon et al. 2005), LIT¬ 
TLE THINGS (Oh et al. 2015), THINGS dwarves (Oh et al. 
201 la, b), THINGS spirals (Oh et al. 2015), Sculptor, For¬ 
nax and the Milky Way. The comparison of the theoretical 
predictions with the complete set of data shows that the two 
models perform similarly, while when we restrict the anal¬ 
ysis to a smaller subsample of higher quality, we show that 
the method presented in this paper (baryonic clumps-DM 
interaction) performs better than the one based on super¬ 
nova feedback. We also show that, contrarily to the first 
model prediction, dSphs of small mass could have cored 
profiles. This means that observations of cored inner pro¬ 
files in dSphs having a stellar mass < 10 6 Mq not necessar¬ 
ily imply problems for the ACDM model. 

Keywords cosmology: theory - large scale structure of uni¬ 
verse - galaxies:formation 
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1 Introduction 

The ACDM model is a highly successful paradigm at large 
scales (Del Popolo 2007; Komatsu et al. 2011; Del Popolo 
2013; Hinshaw et al. 2013; Del Popolo 2014a; Planck Collaboration XVI§ 
2014), but it shows some drawbacks at smaller scales (galac¬ 
tic, and centre of galaxy clusters scales (Del Popolo and Gambera§ 

2000; Del Popolo 2002; Del Popolo and Cardone 2012; 

Newman et al. 2013a, b; Del Popolo 2014a)) 

Of the main problems of the ACDM paradigm the most 
’’stubborn” seems to be the so called Cusp/Core problem 
(Moore 1994; Flores and Primack 1994) 2 dealing with a 
discrepancy between the flat density profiles observed in 
LSBs and dwarf galaxies, and the cuspy density profile ob¬ 
tained in N-body simulations, e.g. the Navarro-Frenk-White 
(NFW) profile (Navarro et al. 1996b, 1997, 2010). The 
NFW profile predicts an inner profile going as p ot r a , with 
a = — 1. An even steeper profile predicted by Moore et al. 

(1998) and Fukushige and Makino (2001) gives p oc r a , 
with a = —1.5, while other authors found that the inner 
slope is dependent on the object considered, and/or its mass 
(Jing and Suto 2000; Ricotti 2003; Ricotti and Wilkinson 
2004; Ricotti et al. 2007; Del Popolo 2010; Cardone et al. 

201 lb; Del Popolo 2011; Del Popolo et al. 2013d; Di Cintio et al.| 

2014). More recent N-body dissipationless simulations 
seem to agree on the fact that a profile flattening towards 


1 For precision sake, the ACDM paradigm suffers of other problems even at 
cosmological scales (e.g., the cosmological constant problem (Weinberg 
1989; Astashenok and del Popolo 2012), and the ’’cosmic coincidence 
problem”). 

-The other most often mentioned problems of the ACDM paradigm, are 
a) the discrepancy between the number of subhaloes that N-body sim¬ 
ulations predict (e.g. Moore et al. 1999) and observations; b) the Too- 
Big-To-Fail (TBTF) problem. In this last problem simulated haloes have 
too many, too dense and massive subhalos in comparison to observations 
(Boylan-Kolchin et al. 2011, 2012). Unified solutions have been proposed 
to the quoted problems, based on the action of baryons located in the inner 
parts of the haloes (Zolotov et al. 2012; Del Popolo et al. 2014). 
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the centre, to a minimum value of ~ —0.8 (Stadel et al. 
2009), namely the Einasto profile, seems to be a better ht to 
simulations (Gao et al. 2008). 

The problem is that the smallest value predicted by dissi¬ 
pationless N-body simulations is larger than the values ob¬ 
tained by observations (Burkert 1995; de Blok et al. 2003; 
Swaters et al. 2003; Kuzio de Naray and Kaufmann 2011; 
Oh et al. 2011 a, b), in SPH simulations (Governato et al. 
2010, 2012), or in semi-analytical models (Del Popolo 
2009; Cardone and Del Popolo 2012; Del Popolo 2012a,b; 


The previous discussion highlights the fact that the de¬ 
termination of the inner slope of galaxies, even dwarves, is 
not at all an easy task. The result from the previous studies 
and several others is that exists a range of profiles, and even 
with the improvements of nowadays kinematic maps there 
is no agreement on the exact dark matter slopes distribution 
(Simon et al. 2005; Oh et al. 201 lb; Adams et al. 2014). 

The situation is even more clear going to larger masses 
(e.g., spiral galaxies) dominated by stars 3 , and especially 
to smaller masses (e.g. dwarf spheroidals (dSphs)) where 


Del Popolo and Hiotelis 2014). Recently, Polisensky and Ricotti|iases that enter in the system modelling (Battaglia et al. 
(2015) found that haloes do not show universal density pro- 2013) lead to opposite results. 


files, rather their shape is determined by the initial linear 
power spectrum of density perturbations. In addition, the au¬ 
thors found that profiles depend on the halo mass, in agree¬ 
ment with previous works on the subject and that warm dark 
matter (WDM) halos develop a core, but this is not signifi¬ 
cant enough to explain observations. 

The cusp/core problem has been also noticed at galaxy 
clusters scales. Kinematics and lensing constraints in cD 
galaxies (BCG) located in the centre of relaxed clusters, 
showed that the clusters DM profiles is flatter than a NFW 
profile, but the total mass profile is in agreement with the 
NFW predictions (Sand et al. 2002, 2004; Newman et al. 
201 3a, b). 

Dwarf galaxies are dark matter (DM) dominated, and 
have a low baryon fraction (de Blok and McGaugh 1997). 
They have been widely used because of their simple dynam¬ 
ical structure, at least disk galaxies without bulges. In the 
case of high-surface brightness objects (larger objects), it 
is more complicated to determine the inner density struc¬ 
ture. So, the previous statement on the cored nature of the 
inner density profile of all galaxies, is not at all obvious. 
While according to Spano et al. (2008) high-surface bright¬ 
ness galaxies are cored, other authors (e.g., Simon et al. 
2005; deBloketal. 2008; Del Popolo and Cardone 2012; 
Del Popolo et al. 2013d; Martinsson et al. 2013) conclude 
differently. The THINGS sample shows a tendency to 
have profiles better described by isothermal (ISO) profiles 
for low luminosity galaxies, Mb > —19 and the pro¬ 
files are equally well described by cuspy or cored profiles 
for Mb < —19. However, even dwarfs do not always 
have flat slopes, as shown by Simon et al. (2005). In the 
case of NGC 2976, 4605, 5949, 5693, 6689, the authors 
showed that the profiles range from 0 (NGC2976) to -1.28 
(NGC5963). Different results have been obtained even us¬ 
ing similar techniques for the same object. For example, in 
the case of NGC2976 the dark matter profile slope is brack¬ 
eted by —0.17 < a < —0.01, according to Simon etal. 
(2003), while a = —0.90 ± 0.15 for Adams et al (2012), 
a = —0.53 ± 0.14 according to Adams et al. (2014), con¬ 
sidering stars as tracers, or a = —0.30 ± 0.18 (Adams et al. 
2014), considering gas as tracer. 


Several techniques have been used. The spherical 
Jeans equation gives results highly dependent on the as¬ 
sumptions, since mass and anisotropy of the stellar orbits 
are degenerate in the quoted model (Evans et al. 2009). 
Maximum likelihood in parameter space in Jeans mod¬ 
elling (Wolf and Bullock 2012; Hayashi and Chiba 2012; 
Richardson and Fairbairn2013) has similar problems. SchwarzschikJ 
modelling has been used for (e.g.) Sculptor and Fornax find¬ 
ing cored profiles (Jardel and Gebhardt 2012; Breddels et al. 

2013; Jardel and Gebhardt 2013; Jardel et al. 2013). Meth¬ 
ods based on multiple stellar populations concluded that 
Fornax (slope measured at ~ 1 kpc) and Sculptor (slope 
measured at ~ 500 pc) have a cored profile (Battaglia et al. 

2008; Walker and Penarrubia201 1; Agnello and Evans 2012;| 
Amorisco and Evans 2012). However, a cusp is found in 
Draco using a Schwarzschild model (Jardel et al. 2013). The 
previous results show that in reality there is no accepted con¬ 
clusion on the inner structures of dSphs. 

On the other side, a clear determination of the cored 
or cuspy structure of dSphs is very important because the 
smaller is the mass of an object the more probable is that 
its inner profile is similar to that of dissipationless N-body 
simulations predictions, namely cuspy. 

Concerning how the Cusp/Core problem could be solved, 
there are at least two different approaches: 

a) cosmological solutions, based on a different spectrum at 
small scales (e.g. Zentner and Bullock 2003), different na¬ 
ture of the dark matter particles (Colin et al. 2000; Goodman 
2000; Hu et al. 2000; Kaplinghat et al. 2000; Peebles 2000; 
Sommer-Larsen and Dolgov 2001), or modified gravity the¬ 
ories, e.g., f(R) (Buchdahl 1970; Starobinsky 1980), f(T) 

(see Bengochea and Ferraro 2009; Linder 2010; Dent et al. 

201 1; Zheng and Huang201 1) andMOND (Milgrom 1983a, b).| 

b) Astrophysical solutions. These are based on the idea that 
the dark matter component of a galaxy expands due to a 
’’heating” mechanism with the result that the inner density 
is reduced. 

The two most known astrophysical solutions are a) ’’su¬ 
pernovae feedback flattening” (SNFF) of the cusp ( | 

3 See Section 3 for a wider discussion. 
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1996a; Gelato and Sommer-Larsen 1999; Read and Gilmore 
2005; Mashchenko et al. 2006, 2008; Governato et al. 2010, 
2012 ), 
and 

b) ’’dynamical friction from baryonic clumps” (DFBC) 

(El-Zant et al. 2001, 2004; Ma and Boylan-Kolchin 2004; 


Mb,in is the initial gas content of the protostructure, while 
Mb is the final baryonic content, namely Mb = M gas+s t a rs • 
Another confirmation of the goodness of the DFBC is that 
the inner slope dependence on the ratio M*/Mh a i 0 is similar 
to SPH simulations results, although in the case of clusters 
and for a different mechanism than that of SPH simulations. 


Nipoti et al. 2004; Romano-Dfaz et al. 2008, 2009; Del Popolo| b) In the cluster case, a further correlation between the inner 
2009; Cole et al. 2011; Inoue and Saitoh 201 1; Nipoti and Binneyjpe and the Brightest Cluster Galaxy (BCG) mass, the core 


2015). 

A wider discussion of the two models is given in Sec¬ 
tions 2.2.1 and 2.2.2. Here, we may recall that the SNFF 
model in the Governato et al. (2010) results has been found 
in agreement with the THINGS galaxies (deBloketal. 
2008; Walter et al. 2008) density profiles (Oh et al. 2008, 
201 la,b). Results somehow contradicting the previous ones 
are those of (e.g.) Trujillo-Gomez et al. (2015) who sim¬ 
ulated feedback from supernovae (SN) and radiation pres¬ 
sure from massive stars, both in disk galaxies and dwarfs. 
They found that the second effect is more important than 
supernovae feedback. The ability of the model to solve 
the small scale problems of the ACDM model has also 
been questioned by several authors (Ferrero et al. 2012; 
Penarrubia et al. 2012; Garrison-Kimmel et al. 2013, 2014; 
Papastergis et al. 2015) (see Section 2.2.1). 

The DFBC predicted the correct shape of galaxy density 
profiles (Del Popolo 2009; Del Popolo and Kroupa 2009) 
in agreement with the SPH simulations performed by 
Governato et al. (2010, 2012) and of clusters (Del Popolo 
2012a) by Martizzi et al. (2012) and as well predicted 
correlations among several quantities observed in clus¬ 
ters of galaxies (Del Popolo 2012a), later observed in 
Newman et al. (2013a,b). 

Recently, Di Cintio et al (2014) showed that the inner 
slope depends on the ratio M*/Mh a i 0 in the SNFF scheme. 
A more detailed discussion of the model is performed in 
Section 2.2.1. It predicts a dependence of the inner slope 
from the ratio M*/Mh a io- The profile goes from a cuspy one 
for low values of the quoted ratio to cored ones and again to 
a cuspy one for large spiral galaxies. We want to recall that 
in the DFBC scheme it was found a correlation among the 
inner slope, the halo mass and the angular momentum of the 
structure (Del Popolo 2009, 2010, 2012a,b), in the case of 
dwarf galaxies, normal spirals, and clusters. In clusters, a 
series of other correlation were found (Del Popolo 2012b), 
namely a) a correlation among the inner slope and baryonic 
content to halo mass ratio at z = 0 , Mb/M 500 , and at Zj n iti a i, 
Mb.in/Msoo (see Del Popolo 2012b, their figures 2, 4) 4 , 
and angular momentum 5 . We recall that the baryonic mass 

4 M 500 is the mass in R 500 , the radius enclosing a density 500 times larger 
than the critical one 

5 Since the angular momentum acquired is inversely proportional to the 
peak height (Del Popolo and Gambera 1996; Del Popolo 2009), the col¬ 
lapse of larger mass structures is slowed down 


radius r core and the effective radius R e , and another correla¬ 
tion between the mass inside 100 kpc, which is mainly dark 
matter, and that inside 5 kpc, mainly constituted by baryons 
(Del Popolo 2014b) was found (see also Del Popolo 2012b). 

Moreover, it was shown that the DFBC mechanism is 
able to explain the flattening in dwarves as well in clusters 
(Del Popolo 2009, 2012a, 2014a). 

In the following, we want to compare the predictions 
of the SNFF mechanism with the DFBC. To this aim, we 
will use the results of Di Cintio et al. (2014), and compare 
those results with the observations of Adams et al. (2014), 
THINGS galaxies, and dwarves, FITTFE THINGS, Sculp¬ 
tor, Fornax, and the Milky Way (MW). 

The paper is organized as follows. In Section 2 we de¬ 
scribe the data used and the models that we compare to the 
data. Sections 3 and 4 are devoted to results and conclu¬ 
sions, respectively. Finally in the appendix A we discuss in 
detail the DFBC mechanism. 


2 Models and data 

In the present paper, we will try to discern among the two 
previously quoted astrophysical mechanisms that are able to 
give rise to cored galaxies. Conscious of the limitations in 
the determination of the inner slope of the density profile a , 
we will use the best data nowadays available. 

2.1 Data 

The data that we will use are those based on ’’high reso¬ 
lution integral field spectroscopy” of seven nearby galax¬ 
ies, namely NGC0959, UGC02259, NGC2552, NGC2976, 
NGC5204, NGC5949, UGC11707, obtained by | 

(2014) through measurements of their gas kinematics and 
integrated stellar light. The complete description of the 
sample selection, photometry, integral field spectroscopy, 
kinematic extraction of gas and stars are discussed in detail 
in Adams et al. (2014). The dynamical parameters of the 
galaxies were obtained using Bayesian statistics, and dif¬ 
ferently from other studies (e.g., dwarf galaxies in Oh et al. 
201 la,b) the entire dark matter density profiles were fitted 
with a Burkert profile (Burkert 1995) 


p(r) = 


Pb 


(1 +r/r b )(l + (r/r b ) 2 ) ’ 


( 1 ) 





4 


Table 1 Galaxy sample properties from Adams et al. (2014). The upper (lower) part refers to the gas (star)-tracer data. 


Galaxy 

a 

logi 0 M, 

Wot 

T» 

log 10 M 200 

log 10 L 



/M e 

km/s 


/Me 

/Lq 

NGC 959 

-0.88 ±0.15 

(0.94+°+ 9 ) x 10 9 

59.86^2.42 

1.10 ±0.15 

11.06 ±0.23 

8.93 

UGC 2259 

-0.72 ±0.09 

(0.25 + qq®) x 10 9 

41-077 

1.07 ±0.27 

11.42 ±0.14 

8.36 

NGC 2552 

-0.38 ±0.11 

(i-3±S:i 7 ) x 109 

65.23l|f 

1.01 ±0.19 

11.33 ±0.11 

9.10 

NGC 2976 

-0.30 ±0.18 

(0.79 ±0.21) x 10 9 

57.13 ±lf 4 

0.83 ±0.22 

11.94 ±0.51 

8.98 

NGC 5204 

-0.85 ± 0.06 

(0.25 ± 0.03) x 10 9 

41.45±i;gg 

1.08 ±0.13 

11.36 ±0.16 

8.37 

NGC 5949 

-0.53 ±0.14 

(2.98tg; 99 ) x 10 9 

82.89177 1 

1.16 ±0.34 

11.82 ±0.42 

9.41 

UGC 11707 

-0.41 ± 0.11 

(1.2±g;| 4 ) x 10 9 

64 4+ 3 -^ 

1.11 ±0.23 

11.49 ±0.18 

9.04 

NGC 959 

-0.73 ±0.10 

(0.92 ± 0.23) x 10 a 

59.55+1“ 

1.08 ±0.27 

11.64 ±0.32 

8.93 

UGC 2259 

-0.77 ±0.21 

(0.25 ±0.1) x 10 9 

41.41+54 

1.10 ±0.44 

11.62 ±0.61 

8.36 

NGC 2552 

-0.53 ±0.21 

(i-6±S:? s ) x 10 9 

69.1lll 0 51 5 

1.24 ±0.55 

11.23 ±0.38 

9.10 

NGC 2976 

-0.53 ±0.14 

(0.89to 2o) x 10 9 

58.981+09 

0.93 ±0.21 

11.56 ±0.46 

8.98 

NGC 5204 

-0.77 ±0.19 

(0.30 ±0.1) x 10 9 

43.671+ 

1.30 ±0.42 

11.76 ±0.51 

9.37 

NGC 5949 

-0.72 ± 0.11 

(3.1 ±0.7) x 10 9 

83.681+°* 

1.20 ±0.28 

11.46 ±0.22 

9.41 

UGC 11707 

-0.65 ± 0.26 

(l-2lg-.gr) x 10 9 

63.78l+« 

1.07 ±0.44 

11.13 ±0.37 

9.04 


and a generalized NFW (gNFW) profile (Zhao 1996): 

ScPc 


p{r) = 


(r/r s ) a [l + (r/r s )] 3 - 


where 


4 = 


200 


3 ((c,a, 1) ’ 


( 2 ) 


(3) 


p c is the critical density and c the concentration parameter. 
The function £(c, a, q h) is defined as (Barnabe et al. 2012): 


C(c,a,q h ) = / 

Jo 


r 2—a 


(1 + r) 


o:—3 


o a/ 1 - (! - 9h) T V c2 


dr 


(4) 


where q^ indicates the 3D axial ratio of the profile. For 
spherical symmetry, q^ = 1 as in Eq. 3. 

The complete list of parameters in the study is listed in ta¬ 
ble 4 of Adams et al. (2014), and apart the DM parameters, 
the gas-based parameters contain a stellar anisotropy term 
/3- that is used as a free parameter in the ’’Jeans Anisotropic 
Multi-Gaussian-Expansion” models (Cappellari 2008). The 
parameters were constrained by using gas and stars as trac¬ 
ers. 

In table 1, we report the galaxy name, the slope a, M *, 
Kot, the mass-to-light ratio T*, M 200 and the luminosity L. 
The stellar mass M, is calculated from the luminosity and 
Y*, while the circular velocity at 2.2 disc scale-lengths Wot 
is calculated by means of the stellar mass Tully-Fisher (TF) 
relation (Eq. 4 of Dutton et al. (2010)): 

lo 'i^ = 2 - 143 + 0 - 28 K log To4^Mj) ' <5) 

where V 2.2 = Wpt . We chose this equation because, as we 
will see later, it was used by Di Cintio et al (2014) in con¬ 
verting M, to Wot (the rotational velocity) in their figure 6. 


Another data set comes from the THINGS dwarfs stud¬ 
ied by Oh et al (201 la, b). The detailed description of the 
parameters of interest to our work and the slope a can be 
found in Oh et al (2008, 201 la, b). Summarizing, they used 
high-resolution HI data from the THINGS survey 6 . They 
selected 7 dwarf galaxies with clear rotation pattern from 
THINGS in order to obtain the rotation curves. In order 
to extract the velocity field from the data cube, different 
techniques can be used (Intensity-Weighted Mean (IWM) 
velocity field; peak velocity fields; single, multiple Gaus¬ 
sian or Hermite polynomial fits). In order to take appro¬ 
priately into account multiple velocity components, non¬ 
circular motions, and so on, Oh et al. (2008) introduced a 
new method to extract from the HI data cube the circularly 
rotating components, the so called ’’bulk-motion extraction 
method”. Since a rotation curve incorporates the dynamics 
of gas, stars and dark matter, in order to obtain the dynam¬ 
ics of the dark matter, it is necessary to extract the baryons 
contribution from the total dynamics. The stellar component 
mass models are obtained deriving the galaxies luminosity 
profiles through a tilted ring modelling applied to the SINGS 
3.6 pm images to obtain the surface brightness profiles. The 
luminosity profiles are then converted into mass density pro¬ 
files using a T* empirical relation, obtained from population 
synthesis models. The HI surface density profile is obtained 
from the column density of HI, and the tilted-ring model ap¬ 
plied to the HI maps gives the radial HI distribution. By 
subtracting the baryons dynamics to the total one, it is pos¬ 
sible to obtain the dark matter mass model of the galaxies. 
The halo models used are the NFW and the ISO profiles. 


6 THINGS was a HI survey program undertaken with VLA comprising 34 
nearby galaxies with high spectral (< 5.2 km/s), and spatial (6") reso¬ 
lution. 3.6 /xm data, with 4" resolution, from SINGS (Spitzer Infrared 
Nearby Galaxies Survey) (Kennicutt et al. 2003), were used to constrain 
the stellar component contribution to the total kinematics. 
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The last profile is given by 


p{r) 


P o 

1 + (r/r c ) 2 


( 6 ) 


where r c is the core radius, and po the halo central density. 

Using different prescriptions for T*, one can obtain a 
’’maximum disk”, a ’’minimum disk”, etc., fit to the rotation 
curve. The density profile can be obtained from the Poisson 
equation (de Blok et al. 2001): 


P (R) 


1 

\ y dV 

(v\ 2 ' 

4t tG 

RdR + 

\Rj _ 


(7) 


Finally, the inner slope is obtained determining the posi¬ 
tion where the slope changes most rapidly (break radius). A 
least-squares fit to the inner points to the break radius (usu¬ 
ally 5 points) gives a. The uncertainty is calculated recalcu¬ 
lating the slope excluding the data point at the break radius, 
and including the first point outside the break radius. The 
error Act is defined as the difference among these slopes. 
Note again that the value of a, in this case, was obtained 
only using the inner points and not the entire density profile 
as in Adams et al (2014). 

In table 2, we show the values of a, M*, and V rot for 
ICG2574, NGC2366, Ho I, Ho II, M81 dwB, DD0153, DDO 
154 (see table 1 of Oh et al. 2011a). 

The same technique previously described can be applied 
to dwarf galaxies from ’’Local Irregulars Trace Luminos¬ 
ity Extremes, The HI Nearby Galaxy Survey” (LITTLE 
THINGS). In table 3, we reproduce a , M*, and U rot for a 
smaller sample of the quoted dwarfs, kindly provided by Se- 
Heon Oh (Oh et al. 2015). 

Finally, in table 4, we present M„ from the THINGS spi¬ 
rals, their V TO t obtained with the Dutton et al. (2010) for¬ 
mula, and the a values obtained from the rotation curve of 
the quoted galaxies. These data were kindly provided by 
Se-Heon Oh. Another determination of the slopes of the 
THINGS galaxies was performed by Chemin et al (201 ) 
(see Section 3 for a wider discussion). 

In our analysis we also used the Simon et al. (2005) 
galaxies not re-studied in Adams et al. (2014) (see table 5). 


Table 2 Galaxy sample properties for THINGS dwarfs from 

Oh et al. (2011a). 


Galaxy 

a 

M, 

(M e ) 

Kot 

(km/s) 

IC 2574 

+0.13 ±0.07 

10.38 x 10 s 

77.6 

NGC 2366 

-0.32 ±0.10 

2.58 x 10 s 

57.5 

Holmberg I 

-0.39 ±0.06 

1.25 x 10 s 

38.0 

Holmberg II 

-0.43 ± 0.06 

2.00 x 10 8 

35.5 

M81 dwB 

-0.39 ±0.09 

0.30 x 10 8 

39.8 

DDO 53 

-0.38 ±0.06 

0.18 x 10 8 

32.4 

DDO 154 

-0.29 ±0.15 

0.26 x 10 8 

53.2 


Table 3 Galaxy sample properties for LITTLE THINGS 

(Oh et al. 2015). 


Galaxy 

a 

k>gio M * 

/Mo 

Vrot 

(km/s) 

DDO 210 

-0.70 ±0.04 

5.602 

6.75 

UGC 8508 

-0.38 ±0.16 

6.477 

11.9 

CVnldwA 

0.03 ±0.27 

6.612 

12.98 

DDO 216 

-0.03 ±1.30 

6.934 

15.98 

WLM 

0.02 ±0.02 

7.090 

17.71 

DDO 70 

-0.48 ±0.02 

7.093 

17.72 

IC 1613 

-0.10 ±0.92 

7.288 

20.14 

DDO 126 

-0.39 ±0.05 

7.356 

21.05 

DDO 133 

-0.11 ±0.16 

7.418 

21.9 

DDO 168 

0.62 ±0.36 

7.710 

26.47 

DDO 101 

-1.02 ±0.12 

7.730 

26.81 

HAR036 

-0.50 ±0.02 

7.764 

27.33 

DDO 87 

-0.01 ±0.48 

7.791 

27.87 

DDO 52 

-0.49 ±0.02 

7.857 

29.16 

DDO 50 

0.10 ±0.41 

7.991 

31.7 

NGC 2366 

-0.34 ±0.10 

8.034 

32.64 

IC 10 

-0.25 ±0.32 

8.072 

33.46 

NGC 3738 

-0.44 ± 0.03 

8.096 

34.89 

NGC 1569 

-0.23 ±0.67 

8.316 

39.14 


Table 4 Galaxy sample properties for spiral THINGS (disks) 

(Oh et al. 2015). 


Galaxy 

a 

logi 0 M, 

/Mo 

14 ot 

(km/s) 

NGC7331 

-1.19 

11.26 

263.2 

NGC3031 

-0.80 

10.91 

209.8 

NGC6946 

-0.70 

10.79 

195.4 

NGC3198 

-0.40 

10.49 

159.9 

NGC3521 

-0.10 

11.09 

235.8 

NGC2403 

-0.70 

9.71 

96.5 

NGC7793 

-0.70 

9.44 

81.1 

NGC4736 

-1.0 

10.35 

146.1 

NGC3621 

-0.9 

10.29 

140.5 

NGC2841 

-1.8 

11.13 

241.95 

NGC2903 

-2.0 

10.21 

133.4 
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Simon et a] (2005) performed an analysis of dwarf and LSB 
galaxies based on Ha high-resolution velocity fields for the 
galaxies NGC 5963, NGC 6689, NGC 4605, and NGC 
5949. In the case of NGC 5963 and NGC 4605, CO ve¬ 
locity fields were studied. In order to avoid the usual prob¬ 
lems connected to the long-slit spectroscopy (Simon et al. 
2003; Swaters et al. 2003) in the RCs determination, they 
used two-dimensional velocity fields. Multiple wavelengths 
velocity fields (e.g., CO, and Ha) were obtained in order to 
further reduce systematic errors. Multi-color imaging was 
also used to improve stellar mass-to-light ratio determina¬ 
tion. This improves the step of modelling and removing 
the stellar disk. By means of the photometry they mea¬ 
sured geometric parameters, and used the routine RINGFIT 
(i.e., tilted-ring modelling) to obtain the radial velocity of 
the system and RCs starting from the velocity fields (see 
Simon et al. 2003, 2005, for a detailed analysis). 

Finally, we show the density slope a, the mass of the 
stellar component M* and the maximum velocity V' roT for 
the Milky Way, Fornax and Sculptor in table 6. 

A further comment regarding the data used is at this point 
necessary. When we only know the stellar mass M* of the 
object considered and we need to infer the rotation veloc¬ 
ity Vrot, we used the relation given by Dutton et al. (2010). 
We therefore checked the reliability of this relation with ob¬ 
jects where both the stellar mass and the rotational velocity 
is known. We noticed though that for some points the ve¬ 
locities inferred differ from the true one and therefore few 
points in the a — M* plot are in a different position with 
respect to the same points in the a — V rot plot. Since this is 
not a systematic effect, we assume nevertheless the validity 
of the Dutton et al (2010) relation acknowledging the fact 
that some uncertainties still hold. 

2.2 Models 

2.2.1 Supernovae feedback flattening 

The models we want to compare data with are the two al¬ 
ready discussed mechanisms of cusp flattening, namely the 
’’supernovae feedback flattening” (SNFF) model (Navarro et al 
1996a; Gelato and Sommer-Larsen 1999; Read and Gilmore 
2005; Mashchenko et al. 2006, 2008; Governato et al. 2010, 
2012; Teyssier et al. 2013), and the ’’dynamical friction 
from baryonic clumps” model (DFBC) (El-Zant et al. 2001, 
2004; Ma and Boylan-Kolchin 2004; Nipoti et al. 2004; 


Romano-Dfaz et al. 2008; Del Popolo 2009; Romano-Dfaz et al.( 
2009; Cole et al. 2011; Inoue and Saitoh 2011; Nipoti and Binney( 
2015). 

The importance of baryons in solving the Cusp/Core 
problem was suggested starting from Flores and Primack 
(1994) and stressed in many following works. The first 
mechanism envisaged was connected to supernovae feed¬ 
back. 

Navarro et al. (1996a) showed that the sudden expul¬ 
sion of baryons into the halo in a single event could flat¬ 
ten the profile. However, Gnedin and Zhao (2002) showed 
that a single explosive event has not sufficient energy to 
form a core, while repeated moderate violent explosions 
could reach the goal (however see Garrison-Kimmel et al. 

(2013) for a different point of view). Mashchenko et al. 

(2006, 2008) showed that in primordial galaxies, random 
bulk motions of gas driven by SN explosions could form a 
core, and a similar model by Governato et al. (2010) found 
the same result. Oh et al (201 la, b) compared the aver¬ 
age slope of THINGS dwarves with the simulations by 
Governato et al. (2010), and Governato et al. (2012) made 
a similar comparison for larger objects, and found a corre¬ 
lation among M* and the inner slope for galaxies having 
M* > 10 6 Mq 7 . The Governato’s papers used the code 
GASOLINE (Wadsley et al. 2004), a N-Body+SPH code 
to simulate galaxies. By means of the ’’zoom” technique 
(Katz and White 1993), the resolution for gas particles was 
Af Pjg as = 3 x 1O 3 * M 0 , and M Pi dm = 1.6 x 1O 4 M 0 for 
DM particles, and the softening 86 kpc. The authors per¬ 
formed a run in which stars formed if the hydrogen density 
was > 100/cm 3 (High Threshold (HT) run), and another in 
which stars formed if the hydrogen density was > 0.1/cm 3 
(Low Threshold (LT) run). 

These simulations, similarly to Di Cintio et al (2014), im¬ 
plement SN feedback through the blast wave SN feed¬ 
back (Stinson et al. 2006), and/or early stellar feedback 
(Stinson et al. 2013). Stars with masses larger than 8 M 0 
deposit an energy of 10 51 erg in the interstellar medium 
(ISM). Even metals are allowed to diffuse between the par¬ 
ticles of gas (Shen et al. 2010). The coupling of energy eject 
(from SN to the ISM is obtained using a coupling coeffi¬ 
cient €esf- In the MaGICC simulations (Stinson et al. 2013) 


7 In Galaxies with M* < 10 ®Mq the supemovae feedback mechanism 

was not able to transform cusps into cores. 


Table 5 Galaxy properties from the Simon et al. (2005) sample. 


Galaxy 

a 

M, 

{Mq) 

Wot 

(km/s) 

NGC4605 

-0.78 ±0.04 

2 x 10 a 

74 

NGC5963 

-1.20 ±0.13 

9.3 x 10 9 

114 

NGC6689 

-0.79 ±0.12 

4.5 x 10 9 

94 


Table 6 Galaxy sample properties for Milky Way (Ascasibar et al. 
2006), Fornax and Sculptor (Walker and Penarrubia 201 1). 


Galaxy 

O 

M, 

{Mq) 

Wot 

(km/s) 

Milky Way 

Fornax 

Sculptor 

-1.03 ±0.04 
— 0.39±°;« 

—0.05/q39 

(6.43 ± 0.63) x 10 iu 

(3.12 ±0.35) x 10 7 
(8 ± 0.7) x 10® 

196.5138 

17.8 ±0.7 

i n q+2.2 

1 ' -°-2.0 
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the fiducial e es f = 0.1. In the following (Sect. 2.2.3), we 
will discuss more in detail the coupling parameters (e.g., 

Penarrubia et al. 2012). 

Di Cintio et al. (2014), using the same code (GASO¬ 
LINE), and similar parameters showed that if M*/Mh a io < 
0.01% the stellar feedback energy is not enough to turn 
cusps into cores, and one expects cuspy profiles similar to 
the NFW or more cuspy. Going up with stellar mass (better 
M*/Mhaio) the profiles become less steep, and when the ra¬ 
tio M*/Mhaio — 0.5% one gets the flattest profiles. For a 
larger M*/Mh a io ratio the deepening of the potential well of 
the galaxies, produced by a larger number of stars, opposes 
the SNFF mechanism and galaxies have more cuspy profiles. 
Di Cintio et al. (2014) used the stellar mass Tully-Fisher 
(TF) relation (Eq. 4 of Dutton et al. 2010) to have predic¬ 
tions on the DM inner slope on the galaxies observed rota¬ 
tion velocity. The model predicts that the galaxies in which 
the cored profiles should be more evident are LSB galax¬ 
ies (in agreement with observations, de Blok et al. 2008; 
Oh et al. 2011a), while in small velocity (mass) dSphs the 
profiles tend to be cuspy. The situation is more compli¬ 
cated for disk galaxies of larger mass (e.g.. Milky Way) 
which are baryon dominated and the uncertainties in the 
disc-halo decomposition are larger. For those galaxies it 
is difficult to have a clear cut on the cuspy or cored na¬ 
ture of the density profile. As already discussed, for larger 
mass galaxies ( V ro t > 150 km/s), de Blok et al. (2008) and 
Martinsson et al. (2013) showed that ISO or NFW fit equally 
well the density profiles. Other studies (e.g., Donato et al. 
2004; McGaugh et al. 2007) reached the conclusion that 
cored profiles describe well the density profile. 

2.2.2 Gas clumps merging 


containing DM and diffuse gas. The proto-structure ex¬ 
pands to a maximum radius and then re-collapses, first in 
the DM component that forms the potential well in which 
baryons will fall. Baryons subject to radiative processes 
form clumps, which collapse to the centre of the halo while 
forming stars (De Lucia and Helmi 2008; Lietal. 2010, 
(see Sect. 2.2.2, 2.2.3)). In the collapse phase baryons are 
compressed (adiabatic contraction (Blumenthal et al. 1986; 
Gnedin et al. 2004)), so making more cuspy the DM pro¬ 
file. The clumps collapse to the galactic centre, because of 
dynamical friction (DF) between baryons and DM, transfer¬ 
ring energy and angular momentum to the DM component. 
The cusp is heated, and a core forms, before stars form and 
stellar feedback starts to act expelling a large part of the gas, 
leaving a lower stellar density with respect to the beginning. 
Feedback destroys clumps 10 soon after a small part of their 
mass is transformed into stars. The mass distribution is then 
dominated by DM. 

The model described is in agreement with El-Zant et al. 
(2001, 2004); Ma and Boylan-Kolchin (2004); Nipoti et al. 
(2004); Romano-Dfaz et al. (2008,2009); Cole et al. (201 1); 
Inoue and Saitoh (201 1); Nipoti and Binney (2015). 

It is the only model able to explain the correct dependence 
of the inner slope of the DM profile over 6 order of mag¬ 
nitudes in the halo mass, namely from dwarves to clusters 
(Del Popolo 2009, 2010, 2012a,b, 2014a). Here we want 
to stress another interesting point. As previously reported 
the model showed that the inner slope is mass dependent 
(Del Popolo 2009, 2010, 2012a,b). Moreover in Del Popolo 
(2012b) and Del Popolo (2014b) were found correlations 
among the inner slope and other quantities with the BCG 
mass and radius, in agreement with results by Newman et al. 
(2013a, b). 


As already reported, the other mechanism able to trans¬ 
form cusps into cores is that proposed by El-Zant et al. 
(2001, 2004) 8 , based on merging gas clumps of 1O 5 M 0 
in the case of dwarves, and 1O 8 M 0 in the case of spi¬ 
rals 9 10 . Energy and angular momentum transfer from clumps 
to DM can flatten the profile, and the process is the 
more efficient, the earlier it happened, when halos were 
smaller. The effectiveness of the process has been con¬ 
firmed by several authors (Ma and Boylan-Kolchin 2004; 
Nipoti et al. 2004; Romano-Dfaz et al. 2008; Del Popolo 


2.2.3 Drawbacks and differences between the models 

As already described, the SNFF model is based on the ex¬ 
pulsion of gas due to SN explosions, while the DFBC model 
is based on dynamical friction between baryons and DM. 
The two models work in a different way. In the case of the 
SNFF model we start from gas which forms stars which can 
then explode into SN if they posses the correct mass. So, 
in this model in order to produce the density profile flatten¬ 
ing, we need a longer and more complex series of events to 


2009; Romano-Dfaz et al. 2009; Cole et al. 2011; Inoue and Sait®hB;h the goal (cusps transformed into cores). In the case of 
2011; Del Popolo et al. 2014; Nipoti and Binney 2015). the DFBC, it is enough that big gas clumps are present to 

More in detail, as shown in Del Popolo (2009); Del Popolo el fflajen the density profile. In a few words, the DFBC model 


(2014), initially the proto-structure is in the linear phase. 


8 For precision’s sake, the concept that large clouds could heat stellar sys¬ 
tems was proposed by Spitzer and Schwarzschild (1951). 

9 In the Nipoti and Binney (2015) simulations, the clump mass was 10 5 — 

10 6 M©, and the total mass 10 9 M©. 


is more ergonomic than the SNFF model. 

In the last years, it has been shown that the SNFF model 
has some drawbacks. Penarrubia et al. (2012) calculated the 

10 This is allowed since star formation is a not efficient process. 
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Fig. 2 The a-V TO t relation compared to observations. Data come from Adams et al. (2014), Simon et al (2005) (excluding the 
galaxies studied by Adams et al. (2014)), the THINGS spirals, the THINGS dwarves, the LITTLE THINGS, Fornax, and Sculptor 
(Walker and Penarrubia 2011) and the MW (Ascasibar et al. 2006). The red, black, and blue lines describe the model of Di Cintio et al. 
(2014), while the green and yellow lines the model of the present paper. The only difference between the left and right panel is due to the 
galaxy sample in Adams et al (2014): gas traced, in the left panel; stellar traced, in the right panel. The green (yellow) curves are our 
equivalent for the red (cyan) curve of Di Cintio et al. (2014). 



Fig. 1 The a-V TO t relation. The dashed lines refer to the calcula¬ 
tions of Di Cintio et al. (2014), and the solid lines are those calcu¬ 
lated in this paper. The dashed lines and the solid ones are differ¬ 
ent by a factor h inverse, coming from the Tully-Fisher relation of 
Dutton et al. (2010), and discarded in Di Cintio et al. (2014). The 
black lines refer to the slope calculated at 0.01 < r/R v i r < 0.02, 
the cyan lines that calculated at 1 < r/kpc < 2 and the red lines 
refer to the slope calculated at 3 < r/e < 10. 


energy that SN must inject into the haloes in order to re¬ 
move the cusp. According to Penarrubia et al (2012), in 
order to transform a cuspy profile into a cored one in MW 
dSphs, an energy in the range 10 53 — 10 55 erg is required. 
The average energy released in a SNII explosion is of the 
order of 10 51 erg (Utrobin and Chugai 2011). So the ex¬ 
plosions of hundreds to several thousands of SNs could in 
principle produce this huge amount of energy. However, 
the low star formation efficiency in dSphs, suggested by 
their luminous satellites abundance, implies that the real 
contribution to the energy could be lower than that needed 
to flatten a profile. Moreover, while the solution to the 
cusp/core problem with the SNFF model needs a large num¬ 
ber of SNs, and so a large star formation efficiency (SFE), 
the solution of another small scale problem of the ACDM 
model, namely the TBTF problem, places an opposite de¬ 
mand on the SFE. In order to eliminate such a tension one 
or more of the following issues should be true: a) a cou¬ 
pling of energy coming from SN II to DM of the order 
of 1. Such a value of the ’’energy coupling” (csn) con¬ 
tradicts observations. In order to describe the metallicity- 
luminosity relation in dSphs, Revaz and Jablonka (2012) 
obtained a value of csn — 0.05, while the value of csn 
used in Governato et al. (2010) is 0.40. Even this large 
value is smaller than what needed to eliminate the quoted 
tension, b) Cusp removal at high redshifts (z > 6 from the 
Sculptor and Fornax cored profiles). Namely, star forma¬ 
tion should peak at redshifts unexpectedly high (z > 6). 
This conclusion is at odds with the fact that star formation 
went on for 12-13 Gyr in Fornax (de Boer et al. 2012a), and 
for 6-7 Gyr in Sculptor (de Boeret al. 2012b). Even if the 
star formation peak was at that redshift, the tension moves 
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to haloes embedding less stars, the formation of cores in 
dSphs with M* < 1O 7 M 0 requires csn — 1. c) A top- 
heavy stellar initial mass function, d) Considerable satellite 
disruption (e.g., by tidal torques). This last issue is promis¬ 
ing in the solution of the previously discussed problem as 
shown by Zolotov et al. (2012); Del Popolo et al. (2014); 
Brook and Di Cintio (2014); Del Popolo and Le Delliou (2014) 

However, we should add that Governato et a (2010) had 
to use a very high star formation threshold to obtain their 
results. Sawala et al (2014) showed that the high threshold 
assumption is not necessary to obtain the results claimed by 
Governato et al. (2010). 

One additional argument of discussion about simulations 
and the TBTF problem is given by Polisensky and Ricotti 
(2014). The authors could largely reconcile the discrep¬ 
ancy between simulations and observations by using the lat¬ 
est cosmological parameters provided by the WMAP and 
Planck team. They concluded therefore that the strong 
tension observed in previous works on the subject was 
mainly due to the incorrect assumption of the cosmolog¬ 
ical parameters. Another source of tension comes from 
the assumed mass of the Milky Way: assuming its most 
recent estimates, the need of baryonic physics necessary 
to decrease the density of the most massive satellites in 
the Milky Way becomes much less compelling. Despite 
this improvement, one more aspect lacks a satisfactory ex¬ 
planation, namely the problem of missing bright satellites 
just outside the Milky Way virial radius, as discussed in 
Bovill and Ricotti (2011a), Bovill and Ricotti (2011b) and 
Garrison-Kimmel et al. (2013). 

On a similar line. Garrison-Kimmel et al. (2013) dis¬ 
cuss general problems of the SNFF model and in partic¬ 
ular problems in solving the TBTF problem. Compar¬ 
ing the effects of blow-outs of different masses (1O 7 M 0 , 
10 8 M q and 1O 9 M 0 ) they found that a single blow-out of 
a fixed mass has more effect in changing the structure of a 
dwarf in comparison with several repeated blow-outs whose 
mass sums to the that of the single blow-out 11 , contrar- 
ily to the SNFF model which requires repeated blow-outs 
(see Governato et al. 2010, 2012; Pontzen and Governato 
2012). From the point of view of the mass, for the high 
resolution simulations of Garrison-Kimmel et al. (2013) to 
have subhaloes density in line with that observed in MW 
dSphs, a quantity of mass equal to 1O 9 M 0 should be 
ejected, which is marginally exceeding the baryon con¬ 
tent of the dSphs. From the energy point of view, since 
the average energy emitted by SNs explosions is 10 51 erg, 
to match the density of a 1O 6 M 0 dSph, the energy of 
40000 supernovae is needed with an efficiency of 100%. 
For six of the nine classical dSphs, this quantity exceeded 

11 However, repeated blow-outs remove mass preferentially from the centre. 


the number of Type II SN explosions expected. More¬ 
over, they find that explosions can flatten the inner slope 
to a > —0.5, never producing a real core with constant 
density, as predicted by the last versions of the SNFF 
model (e.g., Pontzen and Governato 2012). Ferreroetal. 
(2012) and Papastergis et al (2015) arrived to similar con¬ 
clusions, namely the high improbability for the SNFF 
to solve the TBTF problem. Finally, as several authors 
noticed (Choietal. 2014; Laporte and Penarrubia 2015; 
Laporte and White 2015; Marinacci et al. 2014), nowadays 
hydrodynamic al simulations have not the required resolu¬ 
tion to follow the feedback processes which should trans¬ 
form the cusp into a core. 

Finally Oman et al. (2015) and Oman et al. (2016) with 
the use of hydrodynamical simulations investigate the ‘cusp 
versus core’ problem and conclude that this is better char¬ 
acterized as an ‘inner mass deficit’ problem rather than as 
a density slope mismatch. Investigating simulated dwarf 
galaxies and comparing their properties with a selected cat¬ 
alogue of observed galaxies, they find several discrepancies. 

The authors conclude that to solve these discrepancies, it is 
necessary either to change the dark matter physics, that the 
mass profiles of the galaxies giving rise to the ‘cusp versus 
core’ problem are incorrect when inferred from kinematic 
data or that simulations fail to reproduce correctly observa¬ 
tions, as pointed by other authors before. 

Concerning the DFBC model, to our knowledge there are 
not studies on its limits. However, since the effectiveness of 
this process depends on the clumps properties and that of 
the halo, one could speculate about the a) the origin of the 
clumps, and b) the time-scales for the density flattening in 
comparison with the life of the clumps. 

Concerning the first issue, some galaxies (’’clump clus¬ 
ters” and ’’chain galaxies”) at high redshift show clumpy 
structures (e.g. van den Bergh et al. 1996; Elmegreen et al. 
2004, 2009; Genzel et al. 201 1). At t < 2 Gyr, when the 
gas is infalling towards the disc, radiative cooling induces a 
self-gravity instability which leads to clumps formation (e.g. 
Noguchi 1998, 1999; Agertzetal. 2009; Bournaud et al. 
2009; Ceverino et al. 2010), and finally to the formation of 
disc galaxies. 

Concerning the second issue, as shown by Nipoti and Binney| 
(2015), the flattening process, ffl at , happens on the dynam¬ 
ical friction scale ff r i c ~ 1.4t cross /ln A, with t CIOSS ~ 21 
Myr, namely the time-scale is very short. 

At the same time, in order for the mechanism to work, the 
gas clumps should live for a time longer than that needed to 
redistribute DM in the halo centre. The evaluation of this 
life-time is not easy, since it is connected to a poorly known 
subject, namely the process of star formation and feedback 
(e.g. Murray et al. 2010; Genel et al. 2012; Hopkins et al. 
2012; Bournaud et al. 2014). The star formation time must 
be larger than max[t dy n, fcooi], where t dyn = l/y/Gp gas 
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Fig. 3 The a vs. M * relation. Symbols are like in figure 2, ex¬ 
cept for the points representing the two samples from Adams et al. 
(2014): in black we show the gas-traced sample while in blue the 
stellar-traced sample. 


is the dynamical time, and f coo i is the cooling rate. As es¬ 
timated by Nipoti and Binney (2015) (see their Sect. 3.3), 
this should be several tens of megayears. After star forma¬ 
tion starts, if the clumps are destroyed by supernovae, we 
have to wait another ~ 10 7 yr before the first star explodes. 

For precision’s sake, Genzel et al. (201 1) observed out¬ 
flows coming out from the stellar clumps parts of the clump 
clusters, and Genel et al. (2012) suggested that the clumps 
dissipate in a few tens of Myrs (Murray 2011), to 100 times 
this value in high-z clumps (compatible with the redshift of 
flattening). Moreover, Inoue and Saitoh (201 1) found that 
clumps are even more long-lived than what found in the pre¬ 
vious papers. 

Observation of the MW favours the long life thesis of 
clumps. Moreover, estimates of molecular clouds lifetime 
is > 10 s yrs (Scoville 2013). According to the previous 
discussion, the clumps should orbit for several dynamical 
times before being dissipated. 

3 Results 



Fig. 4 Inner slope vs. Ri n . The last quantity, the innermost radius 
point, is defined by means of a least-square fit to the inner data 
point in the case of THINGS dwarves and LITTLE THINGS. In 
the Adams et al. (2014) case, Ri n is obtained by adding the see¬ 
ing to the fibre radius in quadrature. The solid and dashed lines 
represents the theoretical prediction of NFW and ISO halos. The 
black symbols refer to the Adams et al. (2014) data, the 7 THINGS 
dwarves are indicated by the red symbols visible in the figure, and 
the LITTLE THINGS are represented by the brown dots. The 
THINGS SPIRALS are shown with the green points. The blue 
symbols are the de Blok et al. (2001) data (open circles) and those 
of de Blok et al. (2003) are shown with triangles. 


As previously discussed, there are at least two astrophysical 
solutions to the cusp/core problem: a) the SNFF mechanism, 
and b) the DFBC mechanism. In this section, we compare 
the data on the slopes of galaxies spanning a range of stellar 
masses M* = 4 x 10 6 Mq — 1.6 x 10 n Mg, with the two 
quoted models. 

The results of the paper are presented in figures 1-5. 

Figure 1 shows the inner slope in terms of the rotation 
velocity of the total mass. The dashed lines reproduce the 
a- l/ r ot relation calculated by Di Cintio et al. (2014) for a 
measured in different ranges. The red line represents the 
range 3 < r/e < 10, the blue line the range 1 < r/kpc < 2, 
and the black one, the range 0.01 < r/R v ir < 0.02, where 
e is the softening length, and R v \ r is the virial radius. In 
the case of the galaxies simulated by Di Cintio et al. (2014), 
the red line corresponds, for the lowest mass halo to a range 
0.23 kpc < r < 0.78 kpc, and for the largest 0.94 kpc < 
r < 3.13 kpc. The black line 0.60 kpc < r < 1.20 kpc, 
and for the largest mass galaxy 1.30 kpc < r < 2.65 kpc. 

The black line is fitted by (see Di Cintio et al. 2014): 


Ctblack — 0.132 log 
with 

v = 0.84 ( 

1 V1O 9 M 0 


n 2.58 


ol.99 


-0.58 


+ 0.06 


( M* 


0.26 


V10 9 M O/ 

while the blue and red lines are respectively fitted by 
Ckblue = 0.167619 - log [(l 0 *+2.t4248)-0-699049 + 

^qX+ 2.14248) t-56202' 


( 8 ) 


(9) 


( 10 ) 
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and 

C*red = 0.230967 — log 


(10 


X+2.20929\ -0.493541 


+ 


(10 


X+2.20929\ 149315 


(ID 


with X being 


X = log 


0.0702 


10 6.7120 M , 


4.6570 


-0.57912 ^ io- 2 9657 M°- 25589 


0.50674 


( 12 ) 


The equations for a re <j and abiue have been kindly pro¬ 
vided by A. Di Cintio. 

The previous equations give the a-dependence on M*. 
In order to get the dependence on V TO t, we use (Eq. 5), the 
stellar mass Tully-Fisher (TF) relation (Eq. 4 of Dutton et al. 
(2010)), similarly to Di Cintio et al. (2014). Note that the 
three lines in figure 6 of Di Cintio et a] (2014) and the three 
lines of this paper are different by a factor h inverse, coming 
from the Tully-Fisher relation of Dutton et al. (2010), and 
discarded in Di Cintio et al. (2014). 

In figure 2 we compared the theoretical relations ob¬ 
tained by Di Cintio et al. (2014) and that obtained by means 
of our model with observational data. As before, the red, 
blue and black lines represent the slope a calculated at dif¬ 
ferent distances from the galaxies centre by Di Cintio et al. 
(2014). The green line represents the same relation calcu¬ 
lated with our model. They were calculated using the model 
in Appendix A, which is based on Del Popolo (2009) us¬ 
ing a different recipe of gas cooling, and taking into account 
star formation, reionization, and supernovae feedback as de¬ 
scribed in the Appendix A (see also Del Popolo and Hiotelis 
2014). The model was used to study the evolution of proto¬ 
structures and formation of structures of different masses, 
with final halo mass, and stellar masses similar to that stud¬ 
ied by Di Cintio et al. (2014). For each galaxy the stellar 
to halo ratio was calculated, and the slope calculated in the 
same radial bins of Di Cintio et al. (2014) (e.g., see their fig¬ 
ure 3). Converting the stellar mass into halo mass by means 
of the Moster relation (Moster et al. 2013), we can calculate 
the a-M, relation. This can be converted into an a-V 2.2 
relation 12 , through the Tully-Fisher relation (Dutton et al. 
2010, Eq. 4). Our green (yellow) line corresponds to the red 
(blue) one in Di Cintio et al. (2014), their figure 6. Note that 
the turn in the a-V I0 t relation at ~ 25 km/s (in the case of 
our model) originates from the fact that in our model the in¬ 
ner slope of the density profile of a structure is inversely pro¬ 
portional to the angular momentum (e.g. Del Popolo 2009): 
the larger the last the flatter is the profile. When we move 
from normal spiral galaxies to the dSphs region the angular 


12 V 2.2 is the rotation velocity for late-type galaxies, at 2.2 disc scale- 
lengths. V2.2 — Vopt for the quoted late-type galaxies. 


momentum of the structure strongly decreases producing the 
steepening of the inner slope. 

These models are compared with a set of data composed 
by the Adams et al. (2014) galaxies, the Simon et al. (2005) 

galaxy sample (excluding those restudied by Adams et al. 
2014), the THINGS dwarves (Oh et al. 201 la, b), the THINGS| 
galaxies (Chemin et al. 2011; Oh et al. 2015), the FITTFE 
THINGS galaxies, provided by Oh, and Fornax, and Sculp¬ 
tor, whose inner slope was calculated by Walker and Penarrubia§ 
(2011) without adopting a DM halo model, and directly 
from stellar spectroscopic data. They measured the quantity 
T = A log Mj A log r, finding a value of T = 2.6llo'37 f° r 
Fornax, and T = 2.951 q Jg for Sculptor. The relation among 
T and a is given by ctdm < 3 — T (Walker and Penarrubia 
201 1). We added also a point coming from MW, obtained 
from the best fit to the 511 keV emission (Ascasibar et al. 
2006). The figure shows that for V ro t > 100 km/s, only 
4 points over 9 are compatible with the theoretical models. 

For smaller values of the velocity, the Adams et al (2014) 
galaxy sample shows a larger slope in comparison with those 
of other samples (as observed by the same authors). The 
majority of the THINGS dwarf galaxies have slope larger 
than the predictions. The FITTFE THINGS data can be 
used to constrain the models only at small velocities (< 40 
km/s). They agree with the theoretical data, similarly to 
Sculptor and Fornax, but they have large errors. The plot 
shows that the main differences among the SNFF and the 
DFBC model are evident in the velocity range 50 — 100 
km/s. The DFBC predicts steeper slopes, with maximum 
differences around Act ~ 0.2. The other difference is evi¬ 
dent at small velocities. While the SNFF predicts slopes that 
steepen to cuspy, since at small velocities corresponding to 
M*/Mh a io < 0.01% the energy from supernovae feedback 
is not enough to flatten the profile, the DFBC predicts flatter 
profiles. The reason is connected to the different mechanism 
the DFBC is based on. As already reported, in dSphs star 
formation efficiency is low. This means that if we have a 
fixed quantity of gas, the clumps that it forms can act di¬ 
rectly on DM, while in the SNFF mechanism the gas must 
be converted into large mass stars (and as we told the effi¬ 
ciency is low) to explode into SNs that will inject energy in 
the DM. 

The ability of the DFBC mechanism to form profiles flat¬ 
ter than the SNFF mechanism is important since it implies 
that dSphs of low mass are not necessarily cuspy, as pre¬ 
dicted by the SNFF model. This means that the observation 
of cored profiles at small stellar masses (e.g., < 10 6 Mq) 
would not imply that the ACDM model has scarce possibili¬ 
ties to be correct, as it is predicted by the SNFF model (e.g., 
Madau et al. 2014). 

In the left panel of figure 2, we plot the same quantities, 
except that the Adams’ data are obtained from the stellar 
traced observations. 
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Figure 3 shows a in terms of the stellar mass M*. The 
two figures give similar information as the two panels in 
figure 2. However, in this case we did not use Eq. 5 to 
transform M* into V rn t with the result that we have smaller 
uncertainties. Moreover, the plots show more clearly the 
behaviour of the models especially at small stellar masses. 
The plots show that the DFBC model predicts steeper pro¬ 
files in the range 1O 8 M 0 < M* < 1O 1O M 0 with respect 
to the SNFF model. At small stellar masses, M* ~ 1O 4 M 0 
the profile produced by the DFBC is not cuspy, at least it 
does not become steep as a NFW or steeper as happens to 
the SNFF model. Apart the visual difference among the two 
models, we have analysed which one describes better the 
data. We applied a \ 2 to the data and model (see the follow¬ 
ing subsection). 

Since data having not sufficient spatial resolution can 
give rise to larger values of a , in cored systems, it is use¬ 
ful to study the behaviour of the logarithmic slope in terms 
of the data spatial resolution. In figure 4, we plot the log¬ 
arithmic inner slope as a function of the RC’s resolution. 
As Oh et al. (2011a), we plot data coming from THINGS 
dwarves, de Blok et al. (2001) (open circles), de Blok et al. 
(2003) (triangles). We add data from Adams et al (2014) 
and the LITTLE THINGS objects. As the figure shows, at 
low resolution R m ~ 1 the NFW and ISO slopes are al¬ 
most equal (deBloketal. 2001), while at high resolution, 
Ri n < 1, it is possible to distinguish the NFW and ISO 
slopes. The solid and dashed curves represent the theoretical 
NFW and ISO a-R ln relations. While the THINGS dwarves 
and the LITTLE THINGS deviate significantly from the 
NFW predictions, the Adams et al. (2014) data are somehow 
intermediate between the two models, namely intermediate 
between cuspy and cored halos. 

In the following, in order to constrain the two astro- 
physical models able to transform cuspy profiles in cored 
ones, we will first use a) all the galaxies previously dis¬ 
cussed, and then in order to get the best available distri¬ 
bution of slopes, and constraints on the two quoted mech¬ 
anisms, we will b) include only the highest-quality results 
available in literature. We will then include the sets by 
Adams et al. (2014) and Simon et al. (2005) not re-studied 
in Adams et al. (2014) and some of the THINGS galaxies. 

Concerning the case b, we have to choose a sample. From 
the previous discussions, we know that the data sets are 
noteworthy different. The average values of a is different 
for different data sets. The THINGS dwarf galaxies have 

< a >= —0.29 ± 0.07 (Ohetal. 2011b), which is sig¬ 
nificantly shallower than the Adams et al. (2014) data set 
having < a >= —0.67 ± 0.10 (stellar kinematics), and 

< a >= —0.58 ± 0.024 (gas), or those of Simon et al. 
(2005) with < a >= -0.73 ± 0.44. 

The question is what originates those differences. The 
difference among the THINGS dwarves and the Adams et al. 


(2014) result could be due to the fact that the THINGS 
dwarves have smaller stellar masses. However, this point 
seems not to be so important, because the stellar mass of 
the THINGS dwarves is enough to give rise to a core in 
the quoted galaxies, according to the models of SN feed¬ 
back. Oh et al. (201 lb) measured the inner slopes by using 
a power-law fit to the innermost ~ 5 points of the rota¬ 
tion curve (r ~ 1 kpc). These points are obviously the 
most exposed to systematic uncertainties. The slopes mea¬ 
sured in Adams et al (2014) were calculated by fitting the 
entire density profile. Moreover, a large part of the THINGS 
dwarves have evident peculiarities that could bias the slopes 
measurement 13 . 

The THINGS galaxies have some drawbacks. | 

(2008) didn’t derive the slopes of the THINGS galaxies as 
(for the spirals) these would be dominated by the stars, and 
corrections for them would be too uncertain. The only ones 
where they thought the slope would say something about the 
DM were the ones published in Oh et al. (2011a,b) (which 
are dwarves). 

Other authors (Chemin et al. 2011), tried to calculate 
the slopes of the THINGS galaxies. They re-analysed 17 
galaxies, undisturbed and rotationally dominated, in the 
Walter et al. (2008) sample that coincides fundamentally 
with that of de Bloket al. (2008) and found that the mass 
distributions differ from those of de Blok et al (2008). 

Also Ohetal (2015) calculated the THINGS slopes, 
and the results are in many cases in conflict with those of 
Chemin et al. (201 1). In the present paper, we use Oh et al. 
(2015) data for the THINGS galaxies, for compatibility rea¬ 
sons, since we also used Oh’s data for THINGS dwarves and 
LITTLE THINGS. The reason behind it is that it is impor¬ 
tant to use galaxies whose slope is determined in a similar 
way and approximately at the same distance from the centre. 
Hence our choice for using data from Oh et al. (2015) and 
Adams et al. (2014) and neglecting those from Chemin et al. 
( 2011 ). 

This justifies the need to choose accurately the best sam¬ 
ple from the data we have. 

So, we will include the best behaved THINGS dwarves 
(Ho II and DDO 154), the Adams’ galaxies (NGC 959; UGC 
2259; NGC 2552; NGC 5204; UGC 11707) and those of 
Simon et al (2005) not re-studied in Adams. This set has a 
a = —0.673 ± 0.34 almost independent of whether we use 
gaseous or stellar kinematics, and in strict agreement with 
the slopes obtained by Newman et al. (2013a) (see also Del 
Popolo 2014). 

In figure 5, we compare the a-M * relations with the pre¬ 
vious sample. 


13 Holmberg I. Holmberg II, NGC 2366, DD0154, and DDO 53 have not 
trivial kinematic asymmetries, IC 2574, NGC 2366 have not negligible non¬ 
circular motions, and M81dwB has a rotation curve which declines. 
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As reported in the introduction, the cusp/core problem 
has been also noticed at clusters of galaxies scales. It is 
remarkable that the mechanism explaining the quoted ob¬ 
servations, and the correlations found by Newman et al. 
(2013a, b), described in Del Popolo and Cardone (2012) and 
Del Popolo (2014a), is the same explaining the shallow den¬ 
sity profiles over 6 order of magnitudes in the halo mass 
(dwarfs, clusters). In order to explain the galaxy clusters 
density profiles, Martizzi et al. (2012) had to invoke the 
feedback from active galactic nuclei (AGN) on the distri¬ 
bution of gas, finding a cores of similar size in the DM com¬ 
ponent and the stellar central component. 

3.1 Determination of the best model 


As shown in the previous figures, the different models do 
not seem, as one may expect, to reproduce the data perfectly. 
This is understandable since on one side the theoretical mod¬ 
els are very approximate and suffer of many uncertainties 
due to the poor knowledge of the baryon physics involved, 
while on the other side, observational mass and inner slope 
determinations are affected by systematics and difficulties. 

Therefore it is necessary to adopt a statistical approach to 
evaluate which model describes better the observations. To 
do so we use the x' 2 test defined as 


N 


X 


= £ 


(cTth Ct 0 bs,fc) 


fc= 1 


(13) 


where ath and a 0 bs are the theoretical and observed density 
slope respectively and a 2 the error on each measurements. 
This test strictly requires that the error distribution is Gaus¬ 
sian. If this is not the case (being the error bars asymmetric), 
we will use the geometric mean of the two errors. Note also 
that we do not have any free parameter in the models, being 
the models dependent only on the mass of the stellar com¬ 
ponent M* and the numerical constants determined by the 
underlying physical processes considered. 

We will build several likelihood functions: we first con¬ 
sider the whole sample of data points and then the restricted 
subsample described in Section 3. 

In a second step we will consider again all the data points, 
but we will apply a velocity cut-off. In particular, we will se¬ 
lect all the points with maximum velocity V rot smaller than 
50 and 30 km/s. We do not go below this threshold since 
there would be not enough points. With a threshold of 50 
(30) km/s we can use 27 (16) data points out of the 50 used 
in this work. The optimal sample is made of 10 objects. 

We show our results about the reduced chi-square (Xre<j) 
in table 7 using the numerical and theoretical values of the 
inner slope a as in figure 2. In this way we do not need any 
conversion from the velocity V TO t to the stellar mass M*. 
Note that the values of x 2 ( .,i are very high, of the order of 
hundred per degree of freedom, due to the fact that many 


points lie outside the theoretical curves for many as and the 
X 2 statistics is dominated by points with very small error 
bars. We remind the reader that the reduced x 2 is defined as 
Xred = X 2 1 v where v is the number of degrees of freedom 
of the model. 

As it appears evident, all the models perform very sim¬ 
ilarly, with marginal differences not statistically significant 
when we use all the data points, either with gas or stellar 
tracers. The situation is different when we restrict our analy¬ 
sis to the optimal sample. In this case the for the sample 
containing objects whose inner slope was evaluated using 
stellar tracers is two times smaller than the one for the gas 
tracer sample. This shows the higher quality of the obser¬ 
vations. We observe this improvement only for the optimal 
sample since when we apply a cut-off of 50 km/s in Vrot, 
the two reduced x 2 are once again very similar. It is also 
instructive to compare the two different models studied in 
this work. Comparing the models’ prediction for the same 
range of r, we observe that when using the whole sample, 
the two different models behave roughly the same, with no 
particular preference for one or the other. Same conclusion 
when we consider the cut V TO t < 30 km/s, with the model 
of this paper performing a bit better (worse) for the range 
1 < rj kpc < 2 (3 < r/e < 10). When we restrict the 
analysis to the optimal sample, the model of this paper has 
a reduced Xred a Victor of two smaller than the model based 
on the SN feedback, showing therefore a better fit to the data 
points. 

A usual quantity that can be easily evaluated the rela¬ 
tive probability of the models with respect to each other. 
This is defined as the ratio of the likelihood of the models, 
and in terms of the x 2 statistics is expressed as P(a, b) = 
exp [— (x 2 ~ Xb)/2] where a and b represent the two mod¬ 
els, the SNFF and the DFBC, respectively. In other words 
we are asking which model is more likely to fit the data. 
We show our findings in table 8. As expected the values are 
either very small or very big, due to the fact that the x 2 is 
significantly different from one. This is one more confir¬ 
mation that the models fit the current data very poorly, also 
for the optimal sample, despite a significant improvement in 
terms of probability with respect to the full sample. 

We can therefore conclude that while in general the two 
different models studied in this work perform similarly, 
when we restrict to a subsample of data, the model of this 
paper based on dynamical friction performs better than the 
model based on SN feedback even if not at an appreciable 
statistically significant level. 


4 Conclusions 

The aim of this work is to compare the predictions of two 
models thought to be able to solve the cusp/core problem 
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Table 7 Xred values for different samples and models. 

Sample 

0.01 < r/R vir < 0.02 1 < 

r/kpc < 2 

3 < r/e < 10 

Model of this paper 

Model of this paper 





for 1 < r/kpc < 2 

for 3 < r/e < 10 

All (gas) 

139.95 

144.16 

153.30 

128.12 

150.85 

All (stars) 

136.71 

140.88 

149.15 

125.81 

147.61 

Optimal (gas) 

39.67 

41.47 

49.13 

20.78 

30.66 

Optimal (stars) 

21.01 

22.46 

25.39 

8.18 

12.79 

Wot < 50 km/s (gas) 

53.98 

59.85 

77.28 

41.71 

82.59 

Vrot < 50 km/s (stars) 

47.20 

52.96 

68.74 

37.28 

76.30 

Wot < 30 km/s 

70.00 

78.99 

93.76 

57.72 

114.36 


Table 8 Relative probability between the models for different samples. 


Sample 

1 < r/kpc < 2 

3 < r/e < 10 

All (gas) 

1.83 x 10~ 1Y1 

7.97 x 10~ 2Y 

All (stars) 

5.37 x 10 -161 

4.15 x 10~ 17 

Optimal (gas) 

3.63 x 10" 41 

8.19 x 10~ 37 

Optimal (stars) 

1.26 x 10 -28 

2.33 x 10~ 25 

Vrot < 50 km/s (gas) 

3.67 x 10~ iua 

9.39 x 10 2y 

Vrot < 50 km/s (stars) 

2.89 x 10 -89 

4.89 x 10 42 

Wot < 30 km/s 

5.22 x 10 yu 

1.27 x 10 OY 



Fig. 5 Like figure 3, but now the data are only: Adams et al. 
(2014) galaxy sample, Simon et al. (2005) galaxies (excluding 
those studied by Adams et al. (2014)), and the best behaved 
THINGS dwarves (Ho II and DDO 154). 


(Moore 1994; Flores and Primack 1994). The models taken 
into account are the ones based on the SN feedback (see 
Sect. 2.2.1) and on the baryon clumps (see Sect. 2.2.2). 

The first one assumes as main driver of the flattening the 
internal density profile, the action of the SN explosions on 
the surrounding medium causing the expulsion of gas, while 
the second one is based on the idea that energy and angular 
momentum transfer from baryon clumps to the dark matter 
component can flatten the profile and transform a cusp into 
a core. 

We compare the models with a sample of observational 
data built from different sets (see Sect. 2.1 for a detailed 
description). Here we limit ourselves to a short description. 

One of the catalogues adopted is based on the sample of 
seven nearby galaxies by Adams et al. (2014) using both gas 
and stars as tracers, together with three additional objects 
studied by Simon et al (2005). A second data set is based 
on seven THINGS dwarves studied by Oh et al. (201 la,b) 
and a third sample is based on the twenty-one objects of 
the LITTLE THINGS set, kindly provided by Se-Heon Oh 
(Ohetal. 2015). A fourth set, also provided by Se-Heon 
Oh (Oh et al. 2015) includes the THINGS spirals and finally, 
the last set is built with data for the Milky Way, Fornax and 
Sculptor galaxies. 

Analysing figures 3 and 5 (the first with the full set of data 
points, the second with a restricted subsample made of the 
Adams et al. (2014) and Simon et al. (2005) galaxies, and 
the two THINGS dwarves Ho II and DDO 154) and com¬ 
paring the models with a reduced chi-square Xred analysis, 
we see that the whole data set does not favour any of the two 
models in particular, while a restricted analysis to the men- 
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tioned higher quality subsample shows a clear preference for 
the model based on baryon clumps. 

We also showed that dSphs of small mass (M* < 

1O 6 M 0 ) can have cored profiles. For M* < 10 5 Mq the 
slope is ~ —0.6. This is the main difference between the 
SNFF model and the DFBC of this paper. An important 
consequence is that finding a dSphs having M* < 10 6 Mq 
with an inner profile not cuspy, as predicted by the SNFF 
model, is not death hit for the ACDM model. 

The last point clearly shows that the determination of the 
inner structure of dwarf galaxies is of fundamental impor¬ 
tance to determine the nature of the DM. However, as al¬ 
ready reported the inner structure of DM haloes of dSphs 
is still debated, since it is difficult to distinguish cuspy and 
cored profiles (e.g., Strigari et al. 2014). Can this problem 
been solved in the near future? Some authors hinted to 
this possibility. Future observations from the Subaru Hyper- 
Supreme-Camera (Takada 201 0) or GAIA (de Bruijne 201 2) 
have been indicated as a possible way out from the puzzle. 
In reality, even with those observations the problem will not 
be solved except for some larger dwarf galaxies (e.g. Sagit¬ 
tarius, see Richardson et al. 2014). In fact, as previously 
discussed, one method often used to study the density pro¬ 
files is based on the Jeans equations. The method has a 
drawback, a degeneracy between the density profile and the 
anisotropy parameter /3. The direct determination of this pa¬ 
rameter is not possible having just data on the 2D projection 
of stars radius and from the line of sight component of stars 
velocity. Several improvements to the previous case have 
been proposed (see Battaglia et al. 2013; Richardson et al. 
2014). Better results could be obtained by means of the 2+1 
data sets (meaning that we know two of the three position 
coordinates and one of the three velocity coordinates) (see 
Battaglia et al. 2013). Information on the proper motions 
allows the determination of the density slope at half-light 
radius (Strigari et al. 2007). However, the proper motions 
of dwarf galaxy stars is challenging to determine even with 
GAIA (Richardson et al. 2014), with a maximum astromet¬ 
ric accuracy of 7 //as at magnitude V = 10, while it could 
have been possible with an error of ±0.2, determining the 
proper motions of just 200 stars, with the SIM mission 14 , 
which should have had a higher astrometric accuracy than 
GAIA’s (Strigari et al. 2007). 
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14 The Space Interferometry Mission, or SIM was a planned space telescope 
developed by USA. SIM was postponed several times and finally cancelled 
in 2010. 




16 


A Model 

The model used in the present paper was introduced in Del Popolo (2009); Del Popolo and Hiotelis (2014), and then 
applied in several other papers to study the universality of the density profiles (Del Popolo 2010, 201 1), the density profiles 
in galaxies (Del Popolo 2012a, 2014a) and clusters (Del Popolo 2012b, 2014a), and the inner surface-density of galaxies 

(Del Popolo etal. 2013d). 

The model is a semi-analytical model (SAM) that includes an improved secondary infall model (SIM) (e.g., Gunn and Gott 
1972; Hoffman and Shaham 1985; Del Popolo and Gambera 1997; Ascasibar et al. 2004; Williams et al. 2004; Hiotelis and Del Popolo§ 
2006, 2013; Cardone et al. 2011a; Del Popolo et al. 2013b, c). Differently from previous SIMs, the model considers the ef¬ 
fects of non-radial collapse originated by random angular momentum 15 , adiabatic contraction of DM that baryons produce 
(Blumenthal et al. 1986; Gnedin etal. 2004; Gustafsson et al. 2006), dynamical friction effects (e.g., El-Zant et al. 2001, 

2004), random and ordered angular momentum. The model takes into account reionization, cooling, star formation, and the 
supernova feedback (see the following). 

The model follows the evolution of a perturbation starting from the linear phase, expanding with the Hubble flow till the 
phase of maximum expansion (turn-around). In the following phases of collapse and ’’shell-crossing”, it is assumed that the 
central potential varies adiabatically (Gunn 1977; Fillmore and Goldreich 1984). The final profile is given by 


p{x) 


jOtapE m ) r dlnf(xj) ~\ 1 

f{x i) 3 [ d\ng(x { )_ 


(Al) 


where we indicated the initial radius with Xi, the collapse factor with f(x i) = x/x m (xi), and the turn-around radius with 
x m {%\)* given by 


= g{x i) = Xi 


1 + <5i 


Si - (Hr 1 - l) 


(A2) 


In the previous equation, fi; is the density parameter, and 5-, the average overdensity in a given shell. Our model contains 
DM and baryons. Initially, baryons are in the gas phase. Baryon fraction is set equal to the ’’universal baryon fraction” 
/b = 0.17 ± 0.01 (Komatsu et al. 2009) (0.167 in Komatsu et al. 201 1). The baryonic fraction is obtained from the star 
formation processes described in the following. 

In the model, the ’’ordered angular momentum” h (coming from tidal torques of large scale structures on those on the 
smaller scales) is obtained through the tidal torque theory (TTT) (Hoyle 1953; Peebles 1969; White 1984; Ryden 1988; 
Eisenstein and Loeb 1995). The ’’random angular momentum” j is expressed in terms of the eccentricity e = 
(Avila-Reese et al. 1998), where r max is the apocentric radius, and r mi]l the pericentric radius. A correction on the eccen¬ 
tricity is made, according to simulations of Ascasibar et al. (2004), which consider the effects of the dynamical state of the 
system on eccentricity 


^(t'max) — 0-8 



(A3) 


for r max O.lfta- 

The steepening of the profile produced by the adiabatic compression was obtained following Gnedin et al. (2004), and 
calculated using iterative techniques Spedicato et al. (2003), while the effects of dynamical friction were obtained adding 
the dynamical friction force, calculated as described in Appendix A of Del Popolo (2009), to the equation of motions (see 

Del Popolo 2009, Eq. A14). 

Gas cooling, star formation, reionization and supernovae feedback were included as done by De Lucia and Helmi (2008) 
and Li et al (2010) (Sect. 2.2.2 and 2.2.3). 

Reionization, treated as in Li et a (2010), reduces the baryon content, and the baryon fraction changes as 

fb 

/b,halo(*, M vir ) = [1 +0 26M F(z)/M vir ]3 ’ (A4) 

where Mp, is the ’’filtering mass” (see Kravtsov et al. 2004), and as usual the virial mass is indicated as M v - lT . The reioniza¬ 
tion redshift is in the range 11.5-15. 


15 The random angular momentum arises from the system random velocities (Ryden and Gunn 1987). 
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Gas cooling is treated as a classical cooling flow (e.g.. White and Frenk 1991) (see Sect. 2.2.2 of Li et al. 2010). Similar 
results are obtained using the Ryden (1988) treatment. 

The details of star formation are given in De Lucia and Helmi (2008). The treatment of Croton et al. (2006) is used for 
the supernovae feedback. In Di Cintio et al. (2014) it was used the blast wave SN feedback (Stinson et al. 2006). For our 
purposes, the choice of the formalism, even if similar, is not so fundamental. A fundamental difference among our model 
and the SNFF model (e.g., Di Cintio et al. 2014) is that in our case the flattening process happens before star formation, and 
the source of energy is gravitational. Stellar feedback acts when the core is already formed, and disrupts the gas clouds. In 
the SNFF model the flattening process happens after star formation and the source of energy is stellar feedback. 

Concerning these last steps. Gas forms a disc, and the star formation rate is 

ip = 0.03 M si /t dyn , (A5) 

being /. r i y , the disc dynamical time, and M s f the gas mass above a given density threshold, n > 9.3/cm 3 as in Di Cintio et al. 
(2014). The initial mass function (IMF) is a Chabrier one (Chabrier 2003). The amount of stars forming is given by 

AM* = 'ipAt , (A6) 

where At indicates the time-step. 

The quantity of energy injected by SN in the ISM is 

A.Z?sn = 0.5ehaioAM*Vg 2 N , (A7) 

where V S 2 N = //sn L’sn is the energy injected per supernovae and per unit solar mass. The efficiency reheating disc gas 
efficiency produced by energy is fixed at ehaio = 0.35 (Li et al. 2010). psn = 8 x lO -3 /M 0 gives the supernovae number 
per solar mass obtained assuming a Chabrier IMF (Chabrier 2003), and the typical energy released in a SN explosion is 
Esn = 10 51 erg. 

Energy injection in the gas reheats it proportionally to the number of star formed 

AM rehea t = 3.5AM* . (A8) 


The change in thermal energy produced by the reheated gas is given by 

AE hot = 0.5AM reheat K 2 ir , (A9) 


This hot gas will be ejected by the halo if AEg-^ > Alp int . and the quantity is 


AMeject = 


A E sn — A E hot 

0-5KI 


(A10) 


The halo can accrete the ejected material that becomes part of the hot component related to the central galaxy 

(De Lucia et al. 2004; Croton et al. 2006). 

The results of the previous model are in agreement with previous studies on the cusp flattening produced by heating of 
DM by collapsing clumps of baryons (El-Zant et al. 2001, 2004; Romano-Dfaz et al. 2008; Cole et al. 201 1; Inoue and Saitoh 
2011; Nipoti and Binney 2015), and as previously reported predicted the correct shape of galaxy, and clusters density profiles 
together with a series of correlations found in observations like those of Newman et al. (2013a,b). 
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